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Abstract 

In this paper we compute one-loop corrections to masses and couplings in the minimal supersym- 
metric standard model. We present explicit formulae for the complete corrections and a set of compact 
approximations which hold over the unified parameter space associated with radiative electroweak 
symmetry breaking. We illustrate the importance of the corrections and the accuracy of our approxi- 
mations by scanning over the parameter space. We calculate the supersymmetric one-loop corrections 
to the W^-boson mass, the effective weak mixing angle, and the quark and lepton masses, and discuss 
implications for gauge and Yukawa coupling unification. We also compute the one-loop corrections to 
the entire superpartner and Higgs-boson mass spectrum. We find significant corrections over much 
of the parameter space, and illustrate that our approximations are good to 0(1%) for many of the 
superparticle masses. 
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1 Introduction 



Most precision measurements of electroweak parameters agree quite well with the predictions of the 
standard model [Q]. These experiments rule out many possibilities for physics beyond the standard 
model, but they have not touched super symmetry, which evades precision constraints because it decouples 
from standard-model physics if the scale of supersymmetry breaking is more than a few hundred GeV. 
Definitive tests of supersymmetry will probably have to wait for direct searches at future colliders. 

Once supersymmetry is discovered, a host of new questions arise. For example, one would like to test 
the low energy super symmetric relations between the particle masses and couplings by making precision 
measurements of the supersymmetric parameters. One would like to measure the super symmetric masses 
as accurately as possible to shed light on the origin of supersymmetry breaking. Furthermore, one 
would also like to know whether weak-scale supersymmetry sheds any light on physics at even higher 
energies. Indeed, the successful unification of gauge couplings encourages hope that other supersymmetric 
parameters might unify as well. It is important to measure these parameters precisely at low energies so 
that they can be extrapolated with confidence to higher energies. 

It is in this spirit that we present our calculation of one-loop corrections to the minimal supersym- 
metric standard model (MSSM). We define the MSSM to be the minimally supersymmetrized standard 
model, with no right handed neutrinos, and all possible soft-breaking terms. We believe that this mini- 
mal model provides an appropriate framework for analyzing the phenomenology of supersymmetry and 
supersymmetric unification. 

We approach our calculation in the standard fashion associated with precision electroweak measure- 
ments. We take as inputs the electromagnetic coupling at z6ro moniBntiini, ttem — 

1/137.036, the Fermi 

constant, = 1.16639 x 10^^ GeV~^, the Z-boson pole mass, Mz = 91.188 GeV, the strong coupling 
in the MS scheme at the scale Mz, as{Mz) = 0.118, as well as the quark and lepton masses, rrif = 175 
GeV, mfc = 4.9 GeV, and = 1.777 GeV §. 

From these inputs, for any tree-level supersymmetric spectrum, we compute the one-loop VF-boson 
pole mass, M^y, as well as the one-loop values of the effective weak mixing angle, sin e^X, and the DR 

weak mixing angle, s^. We also compute the one-loop corrections to the quark and lepton Yukawa 
couplings, as well as the masses of all the supersymmetric and Higgs particles. 

We work in the DR scheme, and take the tree-level masses to be given in terms of the running DR 
parameters. For each (bosonic) particle, we determine the one-loop pole mass, 

= m2(Q) - 7^e^(M2) , (1) 

where M{Q) is the tree-level DR mass, evaluated at the DR scale Q, and n(p^) is the one-loop self-energy. 
(As usual, Il{p^) depends on Q and on the masses and couplings of the particles in the loop. There is a 
similar expression for the fermion pole mass.) 

In all our computations we include the full self-energies, which contain both logarithmic and finite 
contributions. The logarithmic corrections can be absorbed by changes in the scale Q. Therefore we 
checked our logarithmic results against the one-loop supersymmetric renormalization group equations. 
Since we write our results using Passarino-Veltman functions @ , some of our finite terms are automatically 
correct. As a further check, we verified that our corrections decouple from electroweak observables. 

We present our complete calculations in a series of Appendices. These appendices include the full 
one-loop corrections to the gauge and Yukawa couplings, as well as the complete one-loop corrections to 
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the entire MSSM mass spectrum. While some of these results are not new (the gauge-boson ^, Higgs- 
boson Q and gluino 0, |9|, |lO| self-energies and the gauge-coupling corrections already appear 

in the literature), we include the full set of corrections in order to provide a complete, self-contained and 
more useful reference. 

In Appendix A we write the tree-level masses in terms of the parameters of the MSSM, and in 
Appendix B we define the Passarino-Veltman functions that we use to present our one-loop results. 
In Appendix C we compute the one-loop radiative corrections to the gauge couplings of the MSSM, 
and in Appendix D we write the one-loop corrections to the masses. Where appropriate, we evaluate 
the corrections to the mass matrices to account for full one-loop superparticle mixing. This allows for 
an accurate determination of the masses and mixing through the entire parameter space. Finally, in 
Appendix E we discuss the radiative corrections to the Higgs boson masses. 

The results in the Appendices hold for the MSSM with the most general pattern of (flavor diagonal) 
soft supersymmetry breaking.^ The parameter space is huge because of the large number of operators 
that softly break supersymmetry. Therefore in the body of the paper we illustrate our results in a reduced 
parameter space, obtained by assuming that the soft breaking parameters unify at some high scale. 

The unification assumption is useful because it reduces the size of the parameter space. Moreover, it 
implies certain mass relations that can be tested once supersymmetry is discovered. In addition, for any 
set of parameters, it allows us to determine the unification scale thresholds that are necessary to achieve 
unification. As we will see, the present set of precision measurements is sufficient to begin to constrain 
the physics at the unification scale. 

We implement the unification assumption by solving the two-loop supersymmetric renormalization 
group equations subject to two-sided boundary conditions. At the weak scale, we assume a supersymmet- 
ric spectrum, and for a given value of the ratio of vacuum expectation values, tan /?, we use our one-loop 
corrections to extract the DR couplings gi, g2, 53, At, A;,, and A,- at the scale Mz- We then use the 
two- loop DR renormalization group equations 112] to run these six parameters to the scale MguTj which 
we define to be the scale where gi and g2 meet. 

We require that the soft breaking parameters unify at the scale Mqut- Therefore at Mgut we fix a 
universal scalar mass, Mq, a universal gaugino mass, M1/2, and a universal trilinear scalar coupling, ^o- 
We then run all the soft parameters back down to the scale M| = Mq + AM^^, where we calculate the 
supersymmetric spectrum using the full one-loop threshold corrections that we present in this paper. In 
section 4 we show that this scale is essentially the scale of the squark masses, and that the other scalar 
masses and the Higgsino mass are correlated with it as well. 

We require radiative electroweak symmetry breaking |13], so the CP-odd Higgs mass, uia, and the 
supersymmetric Higgs mass parameter, |//|, are determined in a full one-loop calculation at the scale Mg. 
The sign of fi is left undetermined. We then iterate the entire procedure to determine a self-consistent 
solution. Typically, convergence to an accuracy of better than 10~^ is achieved after four iterations. 

Once we have a consistent solution, we use the results of the Appendices to illustrate the one-loop 
corrections in the reduced parameter space associated with unification. We display results for a randomly 
chosen sample of 4000 points. Our sample is chosen with a logarithmic measure in the range 1 < 
tan/? < 60, 50 < M1/2 < 500 GeV, 10 < Mq < 1000 GeV, and with a linear measure in the range 
—3Mq < ^0 < 3Mg. (The upper limits on Mq and Mx/2 are chosen so that the squark masses are less 
than about 1 TeV. While larger squark masses are certainly possible, they reintroduce the fine tunings 
^ Our results can be readily extended to include inter-generational mixing at the cost of additional mixing matrices. 
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that supersymmetry is designed to avoid.) 

Each of these points corresponds to a (local) minimum of the one-loop scalar potential with the correct 
electroweak symmetry breaking, and each passes a series of phenomenological constraints: We require 



the first- and second-generation squark masses to be larger than 220 GeV |14], the gluino mass to be 
greater than 170 GeV |jl^, the light Higgs mass^to be greater than 60 GeV [|], the slepton masses to be 
greater than 45 GeV and the chargino masses to be greater than 65 GeV |15]. We also require all the 
Yukawa couplings to remain perturbative (A < 3.5) up to the unification scale, and since we assume that 
i?-parity is unbroken, we enforce the cosmological requirement that the lightest super symmetric particle 
be neutral. 

We derive approximations to the radiative corrections that hold with reasonable accuracy over the 
unified parameter space. Where appropriate, we use scatter plots to illustrate the effectiveness of our 
approximations. The approximations consist of two parts. First we identify the most important contribu- 
tions to the one-loop corrections. In most cases these are the loops that involve the strong and/or third 
generation Yukawa couplings. Then we derive approximations to the loop expressions that hold over the 
unified parameter space. 

In the next section, we discuss the radiative corrections to the effective weak mixing angle, sin^ ^eff**) 
and the 1^-boson pole mass, Mw- We illustrate the magnitudes of the different supersymmetric con- 
tributions to these observables. We also discuss the renormalization of the DR weak mixing angle, s^, 
and comment on the way that it affects the gauge thresholds at the unification scale. In section 3 we 
examine radiative corrections to the third generation quark and lepton masses. We illustrate the differ- 
ent contributions and present approximations which hold to a few percent. We also examine Yukawa 
unification and demonstrate the size of the unification-scale Yukawa thresholds. In section 4 we present 
our results for the radiative corrections to the supersymmetric and Higgs boson particle masses. We find 
large corrections to the masses of the light superparticles. We compare our results with those of the 
leading logarithmic approximation and find significant improvements over much of the unified parameter 
space. These corrections are important for unraveling the underlying supersymmetric structure from the 
supersymmetric mass spectrum. 

2 The weak mixing angle and the VK-boson mass 

The calculation of supersymmetric contributions to electroweak observables began in 1984 Q. Since 
then, the precise confrontation of electroweak data with theoretical predictions of the MSSM has been an 



active area of study 11]. Global fits to precision data in the MSSM have been performed by several 
groups [17|. In this section we display our results for two electroweak observables over the parameter 
space associated with radiative electroweak symmetry breaking and universal unification-scale boundary 
conditions. We extract the contributions from the various superpartners, and illustrate the manner in 
which the different contributions decouple from the low energy observables. 



■^The light Higgs boson is similar to that of the standard model in almost all of our parameter space, so we apply the 
standard model bound. 
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Figure 1: Supersymmetric corrections to tlie effective weak mixing angle, sin^ ^1*9* ■ Figure 
(a) shows the corrections from top and bottom squark loops versus the heavy top squark 
mass Wf^; (b) shows the neutralino/chargino contribution versus the light chargino mass; 
(c) shows the slepton correction from all three generations against mgj^ ; and (d) shows the 
complete supersymmetric correction plotted against Tn^+ . 



2.1 Effective weak mixing angle 

We start by considering the effective weak mixing angle, s| = sin^ ^eff**- "^^^ full one-loop calculation is 
presented^ for completeness in Appendix C. The complete result is rather involved; for now we simply 
say that the calculation follows the outline presented above: We take aemj G^, Mz-, as{Mz), and the 
fermion masses as inputs, and compute sf as a function of the supersymmetric masses. 

Because we compute the experimental observable s| in terms of other low-energy observables, its one- 
loop supersymmetric corrections decouple as the supersymmetric masses become larger than Mz- From 
Fig. H we see that sf is especially sensitive to light sleptons, and that the sum of the supersymmetric 
corrections is always negative. We did not plot the Higgs boson and first two generation squark contri- 
butions. They are negligible, less than 1 x 10~^ and 4 x 10~^ in magnitude, respectively. The corrections 
to //-decay and the corrections to the Z-i'^-i^ vertex comprise the nonuniversal corrections to sf. The 
former contributes between —3 and 1.5 x 10~^, the latter between ±1.5 x 10^^, and their sum is in the 
range — 4 to 1 x 10~^. 

With rrit = 175 GeV, we find the standard model value of s| varies between 0.2311 to 0.2315 for Higgs 
masses in the range 60 < m/j < 130 GeV. This is subject to an error of 2.5 x 10~^ from the experimental 
uncertainty in the electromagnetic coupling evaluated at Mz, and to corrections of this same order from 
higher loop effects [18|. Furthermore, increasing rrit by 10 GeV decreases sf by 3.3 x 10~'*. 

These predictions for sf should be compared with the LEP and SLD average^ of 0.23165 ± 0.00024. 
Clearly, the standard model calculation agrees quite well with experiment. The additional contribution 
from supersymmetry can lower the value of sf slightly below 0.2300, or about 6a below the experimental 



■^We do not include the supersymmetric nonuniversal Z-vertex contribution in the Appendix. It is a negligible correction 
in the parameter space we consider. 

^The number quoted here assumes lepton universality. 
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Figure 2: Finite corrections to Mw, in MeV. Figures (a-d) are as in Fig. |l|. 

central value. However, we note that higher-order standard-model corrections, changes in Oem and mt, 
and other precision observables should all be systematically taken into account to delineate which regions 
of parameter space are ruled out by these measurements. We do not attempt such a study here. 

The corrections to s| diminish rapidly as the superpartner masses become heavy. For example, if 
we require m^+, m^+, nih > 90 GeV, we find that s| is shifted by at most —8 x 10^^ relative to the 
standard model value. 
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2.2 VF- boson mass 

We now turn to our second precision electroweak observable, and compute the one-loop correction to the 
W-boson pole mass, Myi/. In Fig. |2| we illustrate some of the finite corrections. The full finite correction 
increases the prediction for My/ by up to 250 MeV. As with s^, the contributions from Higgs bosons and 
the first two generations of squarks are small, less than 12 and 8 MeV, respectively. The nonuniversal 
correction is also small, between —6 and 15 MeV. For large supersymmetric masses, the prediction reduces 
to that of the standard model because of decoupling. 

For mt = 175 GeV we find the standard-model value of Mw in the range 80.39 to 80.43 GeV, with 
an error of ±13 MeV from the experimental uncertainty of the electromagnetic coupling. This is subject 



to additional corrections of the same order from higher-loop effects [18|. If we increase mt by 10 GeV, 
we find a 65 MeV increase in Myr. 

From our calculations we find that the MSSM value for M-^ ranges from 80.39 to 80.64 GeV. These 
numbers can be compared to the current world average, 80.33ib0.15 GeV With the current experi- 
mental error, all of the supersymmetric parameter space lies within 2a of the central value. By the end 
of the decade, the error on M\y is expected to be about 50 MeV. If supersymmetry is not discovered 
by that time, one might think that a much more exacting test could be performed. However, the lim- 
its on the superpartner spectrum will also have increased to the point where the effects on weak-scale 
observables from virtual supersymmetry will be greatly diminished. For example, imposing the limits 
, , m/i > 90 GeV, we find that typically 5Mw < 50 MeV, and at most 5M\y = 100 MeV. 
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Figure 3: (a) The unification-scale correction, £g, necessary to obtain as{Mz) — 0.118, 
plotted versus Mq. (b) The maximum and minimum Sg allowed in minimal SU(5) (top two 
regions) and missing doublet SU(5) (bottom two regions), against Mq. (c) Same as (a) with 
as{Mz) = 0.112. (d) Same as (a) with as{Mz) = 0.124. 



2.3 Gauge coupling unification 



We are now ready to study gauge coupling unification |19]. We start by computing the DR electro- 
magnetic coupling constant, a, and the DR weak mixing angle, s^, as described in Appendix C. The 
DR weak mixing angle is closely related to the effective weak mixing ang le, sj. The main difference is 
that is not an experimental observable, so its radiative corrections involve nondecoupling logarithms 
of supersymmetric masses. However, once we subtract these logarithms we find finite corrections to 
which are quantitatively similar to the corrections to sf shown in Fig. |l[ 

In the context of gauge coupling unification, the finite corrections to are very important [20, [2^. 
They play a significant role in determining the required unification-scale thresholds, which are formally 
of the same order in perturbation theory. As we will see, precision measurements already limit the size 
of these thresholds and place constraints on unified models. 

We determine the unification thresholds as follows: we calculate the full one-loop corrections to a and 
s^, and use them to determine the DR couplings gi and g2- We take as{Mz) = 0.118 from experiment 
[0], and apply the supersymmetric threshold corrections to fix the DR coupling at the scale Mz, 



gl{Mz) as{Mz) 



1 - Aa, 



(2) 



where 



Aa, 



2tt 



3 \Mz 



21n (^] 



-Wi ( "^^'^ 



(3) 



The sum indexed by q runs over the six squark flavors. The constant in (^) transforms the MS coupling 
as{Mz) into the DR coupling 33. We also calculate the one-loop DR Yukawa couplings At(M^), Af,(M^), 
and At-(M^), as described in the next section. 
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We then run the six coupled two-loop renormalization group equations up to the unification scale, 
-^GUT) which we define to be the point where gi and g2 meet. At that scale we define the unification 
threshold, Eg, to be the discrepancy between (73 and the electroweak couplings gi and 52, 

53(Mgut) = 5i(Mgut) (l+Sg) . (4) 



In Fig. |3|(a) we plot the threshold. Eg, versus the squark mass scale, Mg, which we define to be 
M| = + 4M2^2- Fro™ the fi gure we see that for as{Mz) = 0.118, unification requires a negative 
unification-scale threshold correction of between —1% and —3%, depending on the weak-scale supersym- 
metric spectrum. For small Mg, the finite corrections to the gauge couplings are comparable in size to 
the logarithmic corrections; they both decrease £g at small Mg. 

In order for a unified model to be consistent with gauge coupling unification, it must be able to 
accommodate values of ~ —2%. Different unified models give rise to different unification-scale threshold 
corrections. In the minimal |]23[| and missing doublet SU(5) models. Eg depends only on the triplet 



Higgs mass [22|, the same mass that enters the nucleon decay rate formulae. The maximum and minimum 



values of Eg in these two models are shown in Fig. ^(b). The two thin regions correspond to the maximum 
values of Eg, obtained by setting the triplet Higgs mass to 10^^ GeV. The two lar ger regions show the 
minimum values of Eg in each model. These values are found by setting the triplet Higgs mass as small as 
possible, consistent with the bounds from nucleon decay [^]. From the figure we see that it is difficult to 
achieve the necessary thresholds in minimal SU(5), but that missing doublet SU(5) has unification-scale 
thresholds in the right range. 

The unification-scale threshold Eg necessary for unification is directly correlated with as{Mz)- In- 
creasing as{Mz) by 5% increases Eg by about 1%, as expected from the one-loop relation 

^^ = 2^^fe,. (5) 
as{Mz) acuT ^ ^ ^ 

We illustrate this in Figs, ^c-d), where we plot Eg for as{Mz) = 0.112 and 0.124. 

3 Quark and lepton masses 

The full set of radiative corrections to the quark and lepton masses is presented in Appendix D. In this 
section we derive approximations to these formulae, valid for the third generation. 

3.1 Top quark mass 

The top quark mass provides an important input for radiative electroweak symmetry breaking. It receives 
strong and electroweak radiative corrections p6| , 10 1. Our approximation begins by eliminating the small 



electroweak corrections, setting g = g' = and = A;, = 0. We then simplify the resulting expressions 
by setting = because mt is much smaller than a typical squark or gluino mass. 
In this limit, the physical top quark mass is given by 



m = rhtiQ) 



Ami 
1 + 

mt 



(6) 
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where the one-loop correction receives two important contributions. The first is the weh-known gluon 
correction,^ 



mt 



tg 



127r2 



2Bo{mt,mt,0) - Bi{mt,mt,0) 



Q 



31n ^ +5 



mt 



(7) 



Note that this contains a logarithmic and a finite piece; the latter gives a 6.6% contribution]^ The second 
correction comes from the top squark/gluino loops, 



\ rnt J 



9l 



<j Bi{0,mg,mij + Bi{0,mg,mi^) 

Bo (0, m-g , ruf^ ) - (0, , m^^ ] 



sin(20t) 



(8) 



where 6t is the top-squark mixing angle, and 
5o(0,mi,m2) = 



In I — I + 1 + 



m 



M2 



In 



Bi{0,mi,m2) 



In 



M2 



1 



+ 



Inx 



(1-x) 



^(1 — x) Inx 



(9) 



(10) 



with M = max(mi, 7712), m = min(mi, 7712), and x = m\/m\. The full B functions are written in 
Appendix B; the formulae presented here are simplifications that hold when the first argument is zero. 

In Fig. ^ we show the complete correction to the top quark mass as well as the contributions from 
the squark/gluino and electroweak loops. The tree-level mass is defined to be rht{mt). We see that the 
squark/gluino loop contribution can be as large as the gluon contribution for TeV-scale gluino and squark 
masses. The electroweak corrections are small because of cancellations. In the figure we also plot the 
difference between the full correction and our approximation. We see that our approximation is good to 
typically ±1%. 

3.2 Bottom quark mass 

Corrections to the bottom quark mass in the MSSM have received much attention because they can 



contain significant enhanced supersymmetric contributions |28, 29 1. These large contributions play an 
important role in Yukawa coupling unification. Previous studies have included only the enhanced contri- 
butions. In this paper we present our results for the full one-loop correction. Moreover, we systematically 
develop approximations to the supersymmetric corrections. In this way we can see the importance of the 
enhanced contributions relative to the full result. 

The corrections to the bottom quark mass are found as follows. Because the bottom quark is light, 
as (777;,) is large and we must resum the gluon contribution. We start with the bottom-quark pole mass, rrn,. 

^Here and in the following, we implicitly perform DR renormalization, so the 1/e poles are subtracted. 

®We have included the two-loop MS contribution Amt = l.llQjmt p7[, which we assume is close to the DR value. 
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Figure 4: Corrections to tlie top quark mass, versus Mq. Figure (a) shows the full one-loop 
correction. Figure (b) illustrates the correction from the squark/gluino loop; the solid line 
shows the gluon contribution for comparison. Figure (c) shows the electroweak corrections. 
In Fig. (d) we plot the difference between the full one-loop result and the approximation 
given in the text. 



We find the standard-model DR bottom quark mass at the scale rrib using the two-loop QCD correction, 



rhb{mb) 



SM 



nib 



1_ f^V' 
\ rub J 



wherej] 



/ Amb\ 



bg 



5 as{mb) 
3 



+ 12.4 



a 



(11) 



(12) 



V mb ) 3 vr V tt / 

and asijnb) is the five-flavor three- loop running DR coupling. We then evolve this mass to the scale 
Mz using a numerical solution to the two-loop (plus three- loop 0(a^)) standard-model renormalization 
group equations [Q. Taking the bottom quark pole mass rrib = 4.9 GeV, and as(Mz) = 0.118, we find 
the standard-model DR value mb{Mzf^ = 2.92 GeV. 

The final step is to add the one-loop corrections from massive particles, 



mb{Mz) = rhbiMz 



>SM 



1 



/Amb\ 
\ mb ) 



massive 



(13) 



We approximate these corrections as follows. We ignore the small W , Z, Higgs, and neutralino contribu- 
tions. This leaves the squark/gluino and squark/chargino loops. 



\ rrib ) 



massive 



\ mb ) 



+ 



\ rub ) 



(14) 



^We do not know the two- and three-loop corrections to ruf, in the DR scheme. Similarly, we do not know the DR 
three-loop QCD contribution to the running of the strong coupling. In both cases we use the MS values. Alternatively, 
we could have used MS equations to run up to Mz, then convert to DR . The difference between the two approaches, 
Am[,(Mz) < 0.05 GeV, is nearly an order of magnitude smaller than the experimental uncertainty in the bottom quark 
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Figure 5: Corrections to tlie DR bottom quark mass rhb{Mz), plotted versus tan/3. Figure 
(a) shows the full one-loop correction; (b) illustrates the full correction from the bottom 
squark/gluino loops; (c) shows the correction from the top squark/chargino loops. Figure 
(d) plots the difference between the full one-loop result and the approximation given in the 
text. 



We then set = 0. The squark/gluino contribution is again given by (||), with the obvious substitution 
i — > 6. To approximate the squark/chargino contribution, we set g = g' = = Xt = 0, except for terms 
that are enhanced by the Higgsino mass parameter fi or by tan /3. We simplify our expressions by setting 
the chargino masses to M2 and fi, respectively. In this case the squark/chargino loops give rise to the 
fohowing terms, 

+ (/X ^ M2) I , (15) 

where i?o(0, mi, is defined in (^), and q (st) is cos^^ {sin 9t). 

In Fig. |5|we show the corrections to the DR bottom quark mass, mi,{Mz), plotted against tan /3. Figure 
^(a) shows the full one-loop correction, while (b) and (c) illustrate the predominately finite corrections 
from squark/gluino and squark/chargino loops. At large tan/3, the top branches in Figs, (a) and (b) 
correspond to ^ < 0, while for Fig. (c) the bottom branch corresponds to /U < 0. The contributions 
from Figs, (b) and (c) tend to cancel. Because of the cancellations and large corrections, previous 
approximations to the bottom quark mass appearing in the literature can be substantially different from 
the full one- loop result. In Fig. |5|(d) we see that the approximation ( |l4|) typically agrees with the full 
one-loop result to within a few percent. 

3.3 Tau lepton mass 

The corrections to the tau lepton mass are of course much smaller than those of the quarks. After 
resumming the two-loop QED corrections which relate the tau pole mass to the DR running mass at 
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Figure 6: Supersymmetric corrections to the DR tau lepton mass mT-{Mz), plotted against 
tan/3. Figure (a) gives the full one-loop correction; (b) illustrates the difference between the 
full one-loop result and the approximation given in the text. 
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Mz [pC|], we obtain the DR mass mr{Mz) = 1.7463 GeV. We approximate the remaining corrections by 
setting = and keeping only those terms proportional to and enhanced by /u or tan /3. The only 
such terms arise from the chargino loops. They give 
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fiM2 tan /3 



167r2 



(16) 



where -Bo(0, rn-i, m2) is given in (|9|). We illustrate the tau corrections in Fig. ^. The full correction ranges 
from —10% to +6%, while the approximation is good to within a few %. For large tan/3, the top branch 
corresponds to /i > 0. 



3.4 Yukawa coupling unification 

In many supersymmetric unified theories, the DR bottom and tau Yukawa couplings are predicted to 
unify at the scale Mqut To test this hypothesis, one must first extract the running DR Yukawa 

couplings Aft and A,- from the DR bottom and tau masses. Our procedure is as follows. We first use the 
formulae in Appendix D to find DR masses for the bottom and tau. We then use the following relations 
to find the DR Yukawa couplings at the scale Mz, 



1 



ThbiMz) = -^Xb{Mz) v{Mz) COS P{Mz) 

mr{Mz) = ^\r{Mz)v{Mz)cosl3{Mz) 
v2 

We determine the full one-loop DR vev v{Mz) from the relation 



(17) 



M| + 7^e^i^(M|) 



-^[g^Mz) + g'\Mz)y{Mz) 



(18) 



where g and g' are the DR couplings, and is the transverse Z-boson DR self-energy. Alternatively, 
the DR vev v{Mz) can be taken from the following empirical fit. 



v{Mz 



248.6 + 0.9 In ( ^\ 



GeV . 



(19) 
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Figure 7: The unification-scale correction, Sb, tliat is necessary to obtain bottom-tau uni- 
fication with the given values of as{Mz) and mb, plotted versus tan/?. The bottom quark 
pole mass is labeled in GeV. 



This expression gives the correct one- loop vev to an accuracy of better than 1%. 

Once we have the DR Yukawa couplings at Mz, we run them to the unification scale Mqut- At that 
scale we define the unification threshold Eb to be the discrepancy between the couplings, 

Afe(MGUT) = A,(Mgut) (l+Sb) ■ (20) 

Of course, the relation between the bottom quark mass and the DR Yukawa coupling A;, depends strongly 
on the QCD coupling as{Mz)- Therefore we compute Eh assuming that £g has already been chosen so 
that as{Mz) is some fixed value. Such an analysis is illustrated in Fig. 0(a), where we plot £b versus 
tan/? for as{Mz) = 0.118. From the figure we see the well-known feature that Yukawa unification is 
possible with Sb — for small tan/3 (1.2^ tan/? ^1.7) and large tan/3 (15^ tan /3^ 40). In the large 
tan P region we distinguish the two cases, fi < and /x > 0. For /i < we see that Eb is always far from 
zero, in the range —24 to —60%. For /x > we find points which permit bottom-tau unification with 
Eb — 0. However, they are not generic; £b depends sensitively on the parameter space, and varies between 
-20 to +307o. 

The discrepancy Eb is sensitive to the value of as{Mz), as well as the input value for the bottom quark 
pole mass. We illustrate this in Figs. @(b-d), for the {as, rub) values (0.118, 5.2), (0.112, 4.9), and (0.124, 
4.9), with rub in GeV. We note that setting as{Mz) = 0.112 and mb = 5.2 GeV, we find solutions with 
\Eb\ < 0.05 over the whole range 1.3 < tan/3 < 30. 

4 Supersymmetric and Higgs boson masses 

We will begin our analysis of the supersymmetric spectrum by discussing some of its general features. 
We will use some of these features when we derive our approximations for the radiative corrections. We 
will be careful to note when we do, so that one can assess the validity of our approximations in other 
scenarios. 
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Perhaps the most striking consequence of the universal boundary conditions is the fact that they 
produce a low energy spectrum which is tightly correlated with the magnitude of the squark mass scale 
M| = Mq +4M-?,„. In Fig. we show the masses rriu,, mgr, ruA, m-+, and m-+ versus Mg. Prom the 

y ^/^ I— ' Al A.2 

figure we see that the squark masses are nearly equal to Mg, while the other masses are generally within 
a factor of two or three. The exceptions to this degeneracy include the Higgs boson, h, which is always 
light, and the additional possibilities of a light top squark, a light Higgs sector, and/or light gauginos. 
Of course, the gaugino masses are nearly proportional to Mx/2- We find that typically m^o ~ 0.4Mi/2) 
m^o ~ ''^x+ — ^-^^1/21 — 2.4Mi/2) but there are substantial variations from weak-scale threshold 

corrections (and from mixing for the charginos/neutralinos) . We show the ratios and mg/Mi/2 

versus Mq in Pig. |9|. 

In the rest of this section we discuss the one-loop corrections to the masses of the gluino, the charginos 
and neutralinos, the squarks, the sleptons, and the Higgs bosons. In the following expressions for the 
mass corrections we implicitly take the real part of the Passarino-Veltman functions. 



4.1 Gluino mass 



The gluino mass corrections are perhaps the simplest of all the mass renormalizations. They have pre- 
viously been studied in Refs. [Q, |^, |l^]; for completeness we list the corrections in Appendix D. The 
gluino mass corrections arise from gluon/gluino and quark/squark loops. The corrections can be rather 
large, so we include them in a way which automatically incorporates the one-loop renormalization group 
resummation. 
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Figure 9: The ratios (a) to^o/Mi/2, (b) m^o/Mi/2, (c) m~+/Mi/2, and (d) mg/Mi/2 



versus Mq. 



The gluon/gluino loop gives 
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: 2So (Ms , Ms , 0) - Si (Ms , Ms , 0) ; 
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Ml 



+ 5 



(22) 



The quark/squark loop can be simplified by assuming that all quarks have zero mass, and that all squarks 
have a common mass, which we take to be Mg^ , the soft mass of the first generation of left-handed squarks. 
We find 
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^i?i(M3,0,MQ,; 



(23) 



Here 
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(24) 



where M = max{p'^ , m'^) and x 



: /m?. As usual, the full mass renormalization contains logarithmic 



and finite contributions. 

The gluino mass corrections are shown in Fig. In the figure we define the tree-level gluino mass to 
be M3(M3), and we evaluate the one-loop mass at the scale Q = Mq. Because we resum the correction, 
varying the scale from Mq/2 to 2Mq changes the one-loop mass by at most ±1%. 

From the figure we see that the leading logarithmic correction can be as large as 20%, while the finite 
correction ranges from 3 to 10%. The finite contribution is largest in the region where the logarithm is 
largest, so the leading logarithm approximation is nowhere good. On the other hand, the approximation 
we provide typically holds to a few percent. It is off by as much as 6% in the region where the full 
correction is 30%. In this region we expect the two- loop correction to be of order 6%. 
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Figure 10: Corrections to the gluino mass versus nig/Mq. Figure (a) shows the complete 
one-loop corrections; (b) shows the leading logarithmic corrections; (c) shows the finite cor- 
rections; and (d) shows the difference between the full one-loop result and the approximation 
given in the text. 



4.2 Neutralino and Chargino Masses 

The complete set of corrections to the neutralino and chargino masses pl| , |8| is given in Appendix D. 
In this section we present a set of approximations to these corrections. These approximations are more 
involved than those discussed above because there are no color corrections that would dominate the 
results. 

Our approximation is as follows. We start by assuming that |/i| > Mi,M2,Mz- (We find that 
M'^/fj? and M2/ fi^ are less than 0.53; see Fig. |^. We work with the undiagonalized tree-level (chargino 
or neutralino) mass matrix, and correct the diagonal entries only, that is, the parameters Mi, M2, and //. 
This approximation neglects the corrections to the off-diagonal entries of the mass matrices, which leads 
to an error of order (a/47r)M^/^^ in the masses. 

We simplify our expressions by setting all loop masses and external momenta to their diagonal values, 
i.e. we set m^o = Mi, etc. We neglect all Yukawa couplings except and Af,. We also ignore the 
mixings of the charginos and neutralinos in the radiative corrections. This also leads to an error of order 
(a/47r)M^//i^ in our final result. 

We further simplify our expressions by setting all quark masses to zero, and by assuming that all 
squarks are degenerate with mass Mq^ , and that all sleptons are degenerate as well with mass M^-^ . We 
also take the Higgs masses to be ruh = Mz and mn = = ^A- This means that we also neglect 

terms of order (a/47r)M^/m^. 

In this limit, the dominant correction to Mi comes from quark/squark, chargino/charged-Higgs and 
neutralino/neutral-Higgs loops. We find 

(^) = " 1C2 { ^ ' °' ^^^^ ^ + ' °' ^^^^ ^ 

+ ^sin(2/?) ( Bo{Mi,^i,mA) - BoiM,,fJ.,Mz)) 
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+ Si(Mi,^,m^) + Si(Mi,/x,M2)| . (25) 

Since Mz, Mi <C /x, we can simplify this expression by setting Mz = Mi = inside the B functions. 
This gives 

/AMi\ 



V Ml y 327r 



I '^'^OmiMq^ + ^QmxMl^ + Omx^iMz — 2-Bi(0, /i, mA] 



+ ^sin(2/3) (^i?o(0,A^,0) - i3o(0, /x, m^) ) " f } , (26) 

where 5o(0, mi, 777.2) and i?i(0, 7721,7712) were defined in (|9|) and (p^, and 6mi...m2 = ln(M^/Q^) with 
M^ = max(777i, 7712). The form of the finite corrections depends on the assumed hierarchy in the 
low-energy spectrum, but the leading logarithms are always correctly given by the 9 terms. We set the 
first subscript of a term, mi, equal to the external momentum. Note that when the renormalization 
scale equals the external momentum, Q — mi, the theta function reduces to the familiar form, 9m\m2 — 
\n{ml/Q'^)e{ml-Q'^). 

The leading logarithmic corrections are easy to read from Eq. (26).^ Note that the terms proportional 
to sin(2/3) are enhanced by the ratio /x/Mi. These finite terms are completely missed in the run-and-match 
approach because they do not contribute to the beta function. 

In a similar way, we approximate the corrections to M2 from quark/squark and Higgs loops. They 

are 

(^) = - ^{9i?i(M2,0,MQj + 3i?i(M2,0,MiJ 

+ ^sin(2/?) (^fio(M2,/i,mA)-So(M2,^,Mz)) 

+ Bl{M2,^Ji,mA) + Bl{M2,^x,Mz)Y (27) 
Setting Mz = M2 = inside the B functions, we find 

(~M2^) ^ ;^|^ I S^MaAfg^ + 36*^/2^/^^ + OUi^iMz -'^Bi{{),ii,mA) 

+ ^sin(2/3) (^i?o(0,/i,0) - i?o(0,/x,m^)) " t} ' ^^^^ 

There are additional corrections to M2 from gauge boson loops. Because M2 enters both the chargino 
and the neutralino mass matrices, the corrections differ slightly for the two cases. However, to the order 
of interest, it suffices to use the neutralino result, 

VW) = |^|2i?o(M2,M2,Mty)-i?i(M2,M2,M^y) j . (29) 

Because M2 is of order M^y, one must use the full B functions in this expression. Alternatively, one can 
use the following empirical fit which works to better than 1%, 

^2 



/ AMoX^^'"'^" r r3 / Mo \ 

- e{M2 - Mw) 
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0.541n +1.15 

\Mw 



(30) 



The logarithmic part of — 2_Bi(0, /x, tua) in Eq. (Gq) is given by 6'mi, 
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Figure 11: Corrections to the lightest neutrahno mass, as in Fig. 10 



The corrections to fi are obtained in a similar manner. In the hmit g''^ 9^, we find 
3 
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BiifM, M2,mA) + Biifi, M2,Mz) + 2Si(/i, fi,Mz)- 45o(//, /x, Mz) 



(31) 



As above, we set Mz = M2 = inside the B function, in which case ( |3lD reduces to 
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{Xl + A?)Bi(/i, 0, Mq,) + XlBl{^x, 0, Mu,) + XlBiif,, 0, Md, 



^OfiMz - Bi{^,0,mA) + 4 



(32) 



The expression for Bi{p,0,m) is given in Eq. (p4| 



In Fig. 11 we show the corrections to the lightest neutralino mass. In Fig. |ll|(a) we show the full 
correction in percent, with the tree-level mass defined as the eigenvalue of the mass matrix, where the 
running parameters Mi, M2, and /x are evaluated at their own scale. (The tree-level mass matrices also 
contain tan/3 at Mz and the W- and Z-boson pole masses.) The one-loop masses have negligible scale 
dependence. 

As usual, the full corrections are made up of logarithmic and finite pieces. The logarithmic corrections 
are shown in Fig. |TT|(b) and the finite corrections are shown in Fig. pT|(c). Note that the finite corrections 
can be more than half as large as the logarithmic corrections. Indeed, the finite corrections can be larger 
than 5% in the small Mi region, primarily because of the Higgsino-loop term proportional to /i. In 
Fig. |ll|(d) we show the difference between the full one-loop result and our approximation. Here, and in 
the following two figures, the logarithmic corrections [Fig. 0(b)] include an explicit sum over the soft 
squark and slepton masses, while the approximations [Fig. 0(d)] use a single soft squark or slepton mass. 

Figure shows, similarly, the corrections to the lightest chargino mass. Again there is a term 
proportional to fi which dominates the finite corrections when M2 is small. In this region, the finite 
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Figure 12: Corrections to tlie lightest chargino mass, as in Fig. p!o[ 



corrections can be as large as 10%. In Fig. |T2|(d) we show that the difference between our approximate 
correction and the full one-loop mass is less than 2%. These corrections are quantitatively similar to the 
corrections to the second-lightest neutralino mass. 

The corrections to the heavy chargino mass are shown in Fig. These corrections are less than a 
few percent, as are the corrections for the two heaviest neutralino masses. The logarithmic corrections 
are in the range to 2.5%, and the finite corrections are in the range to —3%. Figure 13(d) shows that 



our approximation for the heavy chargino mass generally holds to better than 0.5%. Our approximation 
also works to typically better than 1% for the two heaviest neutralino masses, but it can be off by nearly 



70. 



4.3 Squark masses 

The first two generations of squarks receive QCD [^] and electroweak corrections. However, it is a 
very good approximation to ignore the electroweak graphs, since the dominant corrections come from 
gluon/squark and gluino/quark loops. Neglecting the quark masses, these corrections are as follows. 
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1 -|- 3x -|- (x — 1)^ In |x — 1| — ^2 Inx -|-23;ln 



(33) 



(34) 



and X = m'^/m'^. 

For the case of universal boundary conditions the gluino mass is less than or roughly equal to the 



squark mass, so the correction (|34| ) is essentially finite at Q = rriq. From Fig. 14 we see that it varies 
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Figure 13: Corrections to tlie heaviest chargino mass, as in Fig. nC 
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Figure 14: (a) Full one-loop corrections to the first generation squark mass, niuj^, versus 
the ratio rrig / ■ (b) The difference between the full corrections and the approximation 
in the text, versus nig/muj^. (These are essentially the electroweak corrections.) 



from around 1% for x ^ 1 to between 4 and 5% for x ~ 1. We also see that the electroweak corrections 
are small, less than 0.5%. 

The third generation squark masses receive Yukawa corrections on the order of, and opposite in sign 
to, the QCD corrections. In Fig. 15 we show the full corrections to the third generation heavy squark 
masses. As usual, the tree-level masses are defined in terms of the gauge-boson and quark pole masses, 
as well as the soft masses Mq{Mq), Mu{Mu), and Md{Md)- The tree-level mass matrices also contain 
tw.p{Mz)-, fJ'ifJ'), and Ai{max{\Ai\, Mz)) , where Ai denotes the top or bottom ^-term. (Our convention 
for the third generation squarks is to associate the subscript 1 with the mostly left-handed squark. Since 
the light top squark is predominantly right-handed, its mass is denoted "rn^^.) 

From Fig. ^ we see that the heavy top squark mass receives corrections in the range —5 to 2%, while 
the bottom squark masses receive corrections mostly in the to 3% range. We note that in none of 
these cases does the leading logarithm approximation work well: as is the case for all the squarks and 
sleptons, these corrections are essentially non-logarithmic. (The light top squark mass does receive some 
substantial logarithmic corrections, but they are generally not larger than the finite corrections.) 
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Figure 15: (a) The corrections to the heavy top squark mass versus its mass, (b) The 
difference between the full one-loop heavy top squark mass and the approximation, Eq. (|35|). 
(c) Same as (a), for bi. (d) Same as (a), for 62. 



We will now present our approximation for the top squark mass matrix. We will derive our approx- 
imation for the case of the light top squark, but it also works quite well for the heavy top squark (see 



Fig. 15(b)). The mass of the light top squark receives potentially large additive corrections proportional 
to the the strong coupling and the top and bottom Yukawa couplings. We approximate the corrections to 
m^^ by neglecting g, g' and the Yukawa couplings of the first two generations. We also neglect all quark 
masses except mt, which eliminates all sfermion mixing except for that of the top squarks. 

We neglect the mixing of charginos and neutralinos, so the two heavy neutralinos and the heavy 
chargino all have mass We also make the approximations nih = Mz and mn = mjj+ = rriA- Finally, 
we set p = m the i?-functions if any of the other arguments is much bigger than mj^ • 

This gives the following expressions for the one-loop corrections to the top squark mass matrix, 
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(35) 



The AM^ entries are as follows: 
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+ A(et,/3-|)Bo(0,m,-^,0) + K{et-^,P-^)B^{mi^,mi^,mz] 
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We have defined the two functions 



(36) 



(37) 



(38) 



K[9t,(5) = {2mt cos P cos 6t — sin P — At cos P) sin Ot)'^ 

+ (/i sin P-At cos Pf sin2 6*4 (39) 
n{et,P) = 2mjcos^P sin 26t — 2mt cos P{fi sin P — At cos P) . (40) 

Note that the running mass matrix MtiQ) Eq. ( pSD contains the soft masses Mq and Mjj (as well as 
At, etc.) at some common scale, Q. In the limit A;, — > 0, these expressions are equivalent to the results 
of Ref. [|lO|, with certain external momenta set to zero. 

These approximations depend on the mass of the light top squark. Normally, one would take it to 
be the tree-level mass. For the case at hand, however, the choice is more subtle because for very light 
top squarks, the radiative corrections can be quite large. In fact, the radiative correction can change 
the top squark mass squared from negative to positive. Therefore we shall take m^^ to be the one-loop 
pole mass, which we find by iteration. (We find the one-loop top squark mixing angle by iteration as 
well.) We show the light top squark one-loop pole mass versus the tree-level mass in Fig. [l^(a), where 
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Figure 16: (a) The full one-loop light top squark mass, versus the tree-level mass (in GeV). 
On the X-axis we plot sign(m?_^)|m?^ 1^/^, so a negative tree-level mass corresponds to < 0. 
(b) The difference between the full correction and the approximation in the text, versus the 
one-loop mass, (c) The light top squark mass at one loop, versus Ao/Mg. The solid hue 
corresponds to point (I) in the text, the dashed to point (II). (d) The running and one-loop 
light top squark mass versus the renormalization scale Q, for the choice of parameters (I) 
(solid) and (II) (dashed) in the text, with Ao/Mg — —1.83 and —2.2, respectively. The 
running mass curves each have points where m~ becomes negative. In these cases we plot 
the signed square-root of the mass-squared, as in (a). 



the tree-level mass is the eigenvalue of the mass matrix which contains the running parameters and 
Mq evaluated at their own scale (or Mz, whichever is larger; the tree-level mass matrix also contains the 
top quark and Z-boson pole masses, tan P{Mz), /^(a*), and ^j(max(|At|, Mz))- In Fig. |l^(b) we see that 
our approximation for the light top squark mass holds to within 10 GeV. 

With the present unification assumptions, a top squark with mass less than Mz requires that the 
RR term in the mass matrix be small and that the LR element, proportional to the ^d-term, be large. 
The light top squark mass results from a cancellation between the diagonal and off-diagonal terms, which 
requires a fine tuning. We illustrate this in Fig. p^(c), where we plot the light top squark one- loop mass 
versus Ao/Mg. On the same plot we show the curves corresponding to two choices of parameters, (I) 
tan/3 = 20, Mo = 500 GeV, M1/2 = 100 GeV, and < 0, and (II) tan/? = 5, Mq = 100 GeV, M1/2 = 200 
GeV, and /i > 0. Whether at tree-level or one-loop, the parameter ^0 must be tuned to one part in 75 
to obtain a light top squark mass below 50 GeV. 

We see from Fig. 16(a) that the light top squark mass-squared can be raised from —(100 GeV)^ at 
tree-level to over (100 GeV)^ at one loop. For such large corrections, it is important to keep in mind that 
two-loop effects might be important. The size of these effects can be estimated by the scale dependence of 
the one-loop mass. In Fig. [l^(d) we show the scale dependence at the points (I) and (II), with ^0 = —985 
GeV and ^0 = —907 GeV, respectively. We see that as the renormalization scale increases from 100 to 
1000 GeV, the running masses vary over a wide range. In contrast, the scale dependence of the one- loop 
masses is quite mild. 
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Figure 17: (a) The full set of corrections to the left-handed selectron mass versus its mass, 
(b) The complete corrections to the electron sneutrino mass versus its mass. The complete 
corrections to the predominantly (c) left-handed and (d) right-handed tau slepton masses. 



4.4 Slepton masses 

Corrections to the SU(2) mass sum rules were studied in Ref. |3^. They suggest that the corrections to 
the slepton masses are small. We find that the corrections to the left-handed electron or muon slepton 
masses are typically in the range ±1%, as illustrated in Fig. [T7|(a). The right-handed electron and 
muon slepton mass corrections are larger, but still less than 1.7%. The sneutrino mass corrections are 
essentially identical for all three generations. The full correction is typically in the range ±1%, and 
reaches at most —2.5% at ~ 50 GeV, as shown in Fig. [T7|(b). The (predominantly) left- and right- 
handed tau slepton corrections are similar to the corrections of the first two generation charged slepton 
masses. However, from Figs. |l^(c-d) we see that the scatter plots show less uniformity because of the 
additional Yukawa coupling corrections. We emphasize that in the corrections to the slepton masses, 
the leading logarithmic approximation [^ ] typically gives zero correction. (The mass ruf^ receives some 
logarithmic corrections. However, the finite corrections are typically of the same order as or larger than 
the logarithmic corrections.) 

4.5 Higgs boson masses 

We first discuss the corrections to the heavy Higgs boson masses {niA, mn, mfj+), and then we consider 
the one-loop light Higgs boson mass. The full one-loop corrections to the Higgs boson masses appear in 
Ref. [|. 

As usual, we parametrize all the Higgs boson masses at tree-level in terms of the CP-odd Higgs boson 
mass, ruA, and tan/3. To compute the one-loop mass niA, we first take the soft masses m'j^_^{Q) and 
iQ) outputs from the renormalization group equations, and apply corrections from the electroweak 
symmetry breaking conditions to obtain the DR running parameters m\{Q) and fJ-'^{Q) (see Appendix 
E). We then apply further corrections to obtain the CP-odd Higgs boson pole mass, m^, from the running 
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Figure 18: (a) The one-loop pole mass, toa, versus the running mass at the scale niA (the 
ordinate is the signed square-root of the mass-squared, sign{m\)\m\\^ /^); (b) the corrections 
to the charged Higgs mass, m//+ , versus m//+ . 



mass, TTiAiQ), 

m\ = rh\{Q) - neI[AA{m\) + 4^ + si- , (41) 

^ V2 ^ f 1 

where ti and t2 are the tadpole contributions listed in Appendix E. Note that we treat the Higgs mass 
analogously to the superpartner masses, in that we compare the pole mass with the running mass. 
However, the Higgs mass is different because the tadpole, or effective potential, corrections must be 
added to the "tree- level" running mass to obtain the DR running mass, iJiAiQ)- 

The difference between the running mass and the pole mass can be quite substantial; as for the light 
top squark, the radiative corrections can change a negative mass-squared running mass into a positive 
mass-squared pole mass. In Fig. |l8|(a) we show the one-loop pole mass, m^, evaluated at Q = Mg, versus 
the running mass, m^(myi). 

From Fig. [l8|(a) we see that there are points in parameter space where the running mass-squared is 
— 1 TeV^, while the one- loop mass is over 300 GeV. (If one considers the running mass at the scale Mz, 
the largest corrections are even more extreme.) These large corrections arise from terms enhanced by 



tan (3. In fact, all of the points in Fig. 18(a) with m'^{mA) < occur for tan/3 > 25. The tan/3 enhanced 
contributions come only from the last term in (|4l| ) since 1/vi scales like tan/3 for large tan/3. Therefore 
the tan /3 enhanced corrections are simply 

Q 9 35/a^tan/3 f q q 

- rh^{Q) = — — | A^^t 5o(0, m^-^ , m^-J + Af, A -Bo(0, m^^ , m^J 

+ /M2 (^1 - So(0,/i,M2)| , (42) 

where i?o(0, mi, m2) is given in These terms account for the large corrections seen in Fig. [l^(a). 

We parametrize the other heavy Higgs boson masses in terms of the pole mass, niA- The corrections 
to the heavy Higgs boson masses, mj^+ and rriH, turn out to be quite small. For example, the corrections 
to mjj+ are typically less than 1%, as shown in Fig. ^(b). Corrections to the charged Higgs mass are 
the subject of Ref. [^] . 

The corrections to the light Higgs boson mass, m|, have been studied extensively in the literature 
H, ||,|3§. Here we show our results for the full one-loop pole mass over the parameter space associated 
with radiative electroweak symmetry breaking and universal unification-scale boundary conditions. We 
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Figure 19: Figure (a) shows the full one-loop light Higgs mass versus tan/3. The line 
indicates the tree-level bound nih < Mz\cos2f3\. Figure (b) shows the contribution from 
the gauge/Higgs/gaugino/Higgsino loops; (c) shows the difference between the full one- 
loop result and Dabelstein's approximation; and (d) shows the difference between the full 
one-loop mass evaluated at the scales Mq and Mq/2. 



also show the size of a set of corrections which are often neglected, and the accuracy of Dabelstein's 
approximation p7[ |. 

In Fig. |l^(a) we show the one-loop light Higgs boson mass versus tan (3, as well as the upper limit at 
tree level, ruh = Mz \ cos 2/3|. The upper limit on the one- loop Higgs mass depends sensitively on the top 
quark mass, and somewhat less sensitively on the squark masses. In the parameter space we consider the 
squark masses are less than about 1 TeV. With rrit = 175 GeV, we find nih < 130 GeV. 

In Fig. |l9|(b) we show the gauge/Higgs/gaugino/Higgsino contribution to the Higgs mass, which is 
typically —2 GeV. In Fig. ^9|(c) we show the difference between the one-loop result and Dabelstein's 
approximation, Eq. (4.9) of Ref. which only includes the top sector. We see that this approximation 
is typically 2 to 6 GeV larger than the full one-loop mass. 

In any pole mass, the scale dependence formally cancels. However, at any given order, there are 
usually higher-order corrections which do not cancel. For example, when we vary the scale in the Higgs 
mass calculation, we change the tree-level DR mass. To one-loop order, this variation is canceled by the 
change in the self-energy. However, as the scale varies, the couplings and masses in the self-energies also 
change. For the case of the Higgs mass, the change in the top quark Yukawa coupling gives rise to a 
two-loop 0{Xf) scale dependence in our one-loop results. In Fig. |l^(d) we show the difference between 
the full one- loop result evaluated at the scale Mg and Mq/2. We see that the scale dependence is usually 
small, in the range to —1 GeV. The fact that it is negative follows from the dependence of tan/3 on Q, 
and from the dependence of on tan (3. 
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5 Conclusions 



In this paper we computed one-loop radiative corrections in the minimal supersymmetric standard model. 
We took as inputs the electromagnetic coupling at zero momentum, aem, the Fermi constant, G^, the 
Z-boson pole mass, Mz, the strong coupling in the MS scheme at the scale Mz, as{Mz), and the quark 
and lepton masses. From these we computed the TF-boson mass, M\y, as well as the one-loop value of 
the effective weak mixing angle, sin^ ^eff**' ^ * function of the supersymmetric parameters. 

We studied the size of the corrections in a reduced parameter space associated with the unification 
of the soft breaking parameters and radiative electroweak symmetry breaking. We found that super- 
symmetric radiative corrections can reduce sin"^ 9^^^^ by as much as —1.6 x 10^'^ with respect to the 
standard-model value. Similarly, we found that they can increase Mw by as much as 250 MeV. Because 
of decoupling, the points with the largest deviations are also the points with the lightest superpartners. 
As direct searches increase the limits on the superparticle masses, the size of the supersymmetric radia- 
tive corrections will decrease. Indeed, if superparticles are not discovered at LEP 2, we found that the 
maximum size of the supersymmetric radiative corrections will be reduced by a factor of two. 

The apparent unification of the SU(3), SU(2), and U(l) coupling constants is a major piece of evidence 
in favor of supersymmetry. At next-to-leading order, the weak- and unification-scale threshold corrections 
come into play. The weak-scale thresholds decrease the one-loop weak mixing angle. This leads to an 
increase in the predicted value of the strong coupling, as{Mz)- As we have seen, for squark masses less 
than one TeV, a unification-scale threshold of —1 to —3% is necessary to bring as{Mz) into accord with 
experiment. 

The size of the unification-scale thresholds places an important constraint on unified model building. 
In any unified model, the unification-scale thresholds can be calculated as a function of the grand uni- 
fication parameters. One can see whether the model is consistent with a unification-scale threshold of 
about —2%. In this paper we studied the minimal SU(5) model and the missing doublet SU(5) model. 
We found that the former was not compatible with gauge coupling unification, while the latter was. 

Grand unified theories also predict the unification of certain Yukawa couplings, and in a similar 
fashion, the mismatch of the Yukawa couplings at the unification scale can be used to constrain unification- 
scale physics. To this end it is necessary to extract as precisely as possible the DR Yukawa couplings 
from the fermion pole masses. In this paper we presented full one-loop relations between the two, as well 
as approximations that work at the 0(1%) level. We studied the substantial (up to 50%) tan /3-enhanced 
corrections to the bottom quark mass, as well as the corrections to the top and tau masses, which are of 
order 5 percent. 

Supersymmetry also predicts relations between the masses and couplings of the supersymmetric par- 
ticles. Indeed, if new particles are discovered at future colliders, it will be necessary to check these 
relations to see whether the new particles are in fact supersymmetric partners |38]. The radiative correc- 
tions to the supersymmetric mass spectrum presented in this paper will be an essential element in these 
determinations . 

The corrections to the supersymmetric mass spectrum will be used in (at least) two ways. First, they 
will be used to correct the tree- level mass sum rules [33, 3E] which test supersymmetry at the weak scale. 
Second, they will be needed to extract the underlying soft parameters from the physical observables. The 
soft parameters can then be run to higher scales, to test for unification and possibly to shed light on the 
origin of the supersymmetry breaking. 
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The corrections to the supersymmetric masses in the spin-1/2 sector include potentially 30% correc- 
tions to the gluino mass, as well as 0(5%) corrections to the neutralino and chargino masses. In the 
spin-0 sector, the famous quadratic divergences give rise to large corrections to the scalar masses. These 
corrections can lift the running mass-squared of the light top squark from —(100 GeV)^ to (100 GeV)^. 
Even more dramatically, large tan /? corrections can lift the mass-squared of the CP-odd Higgs boson 
from, e.g., -(1 TeV)^ to (300 GeV)^. 

Radiative corrections also have an important effect on the mass of the lightest Higgs boson, h. In 
the parameter space we consider, they effectively change the sign of the tree-level bound from mh < 
Mz \ cos2/3| to rrih > Mz \ cos2/3|. We found the light Higgs mass was raised to at most 130 GeV. The 
corrections to the rest of the scalar masses are smaller. For example, we found 1 to 5% corrections to the 
first two generation squark masses, and 0(1%) corrections to the slepton masses. 

In the paper we presented approximations to many of the formulae for the supersymmetric mass cor- 
rections. These approximations, often good to better than a couple of percent, provide useful substitutes 
for the full corrections. 
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Appendix A: Tree-level masses 



In this appendix we define the tree-level masses. These tree-level relations also hold for the running 
DR parameters at a common scale, Q. For the most part we follow the conventions of Ref. [^ ]. 

The up- and down-type quark and charged-lepton masses are related to the Yukawa couplings and 
the vev's vi and V2 by 



The ratio of vev's vijvx is denoted tan/3. The tree- level gauge boson masses are 



w 



1 2 / 2 
75 (^1 



^ f 12 , 2\ ( 2 , 2 
j[9 +9 )[Vl+V2 



(A.l) 



(A.2) 



M 



(A.3) 



where g and g' are the SU{2) and U{1) gauge couplings. 

The Lagrangian contains the neutralino mass matrix as —ip^'^Ai^oi'^ + li-C-j where ip'^ = {—ib, 
hi, 112)^ and 

/Ml -Mzcpsw Mzs/3sw \ 

M2 Mzcpcw -Mzspcw 

-Mzcpsw MzCfscw fj, 

\ MzSfssw -MzSfscw fJ^ / 

We use s and c for sine and cosine, so that s/3 = sin/3, C2/3 = cos 2/3, etc. Mi and M2 are the soft 
supersymmetry-breaking bino and wino gaugino masses, /U is the supersymmetric Higgsino mass, and 
sw (cw) is the sine (cosine) of the weak mixing angle. The neutralino masses are found by acting on 
the matrix M^o with a unitary matrix A^, so that N*AA^qN'^ is a diagonal matrix which contains the 
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physical neutralino masses, m^o. In the usual case that one of the eigenvalues of (A. 3) is negative, the 
matrix N is complex even if the elements of J^^o are real. 



The Lagrangian contains the chargino mass matrix as —ip '^Ai^+ip'^ + h.c, where = {—iw~^, ^2)'^ ^ V' 
{—iw~, hi)'^ and 

f M2 V2MwSf3 \ 



M 



1/)+ 



(A.4) 



The chargino masses are found by acting on the matrix A4^+ with a biunitary transformation, so that 
U*A4j^+V^ is a diagonal matrix containing the two chargino mass eigenvalues, m-+. The matrices U and 

V are easily found, as they diagonalize, respectively, the matrices A^^_^A^T^ and A^t^A^^+. 
At tree level the gluino mass, rrig, is given by the soft mass, M3. 
The tree-level squark masses are found by diagonalizing the following mass matrices, 

\ ruu {Au + cot (3) 



nid {Ad + /itan/3) 



m„ {Au + /ucot (3) 
M^ + ml + gu^Mlc2f3 

TTid {Ad + /itan/3) 
Ml + ml + gd^Mlc2p 



(A.5) 
(A.6) 



Here Mq, Mjj, and are the soft supersymmetry-breaking squark masses, and the Af^s are the soft 
supersymmetry-breaking A-terms. The slepton mass matrices are analogous. The soft slepton masses 
are denoted and Me- We have defined the weak neutral-current couplings 
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(A.7) 



The electric charge, hypercharge, and third component of isospin of the sfermions are^ 





UL 






d-R 


9 








2 
3 


2 
3 


1 
3 


1 
3 





-1 


1 


Yf 


1 
3 


4 
3 


1 

3 


2 

3 


-1 


-1 


2 


H 


1 
2 





1 
2 





1 
2 


1 
2 






(A.8) 



A symbol without an L or subscript refers to the L-field (e.g. = 2/3). 
The matrix which diagonalizes a sfermion mass matrix is denoted by 

( Cf Sf\ 



where cj is the cosine of the sfermion mixing angle, cosOf, and Sf the sine. These angles are given by 

tan(2ej = 2m {A cot P) 

^ Ml-M^ + {\-2euS^^)Mlc2p' ^ ^ 

tan(20,) = 2rn,(^, + Mtan^) 

^ Ml^-Ml + {-\-2edslr)Mlc2p ^ ' 



Our convention for the right-handed sfermion fields is the charge-conjugate that of Ref. 
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Since there is no right-handed sneutrino, the slepton mixing angle for u satisfies Cy = 1, and the sneutrino 
mass is m? = + M|c2/3/2. 

Given values for tan/3 and the CP-odd Higgs-boson mass, m^, the other Higgs masses are given, at 
tree level, by 



and 



(A.12) 



(A.13) 



The CP-even gauge eigenstates (si, S2) are rotated by the angle a into the mass eigenstates (i?, K) as 
follows, 

\ h ) \ -Sa Ca J \ S2 J 



At tree level, the angle a is given by 



tan 2a = i tan 2/3 . 

ml -Ml 



(A.15) 



Appendix B: One-loop scalar functions 

The following integrals appear at one loop in a self-energy calculation 

1 

i (27r)" q"^ — m? + ie 
(Pq 1 



Bo{p,mi,m2) = 167r2Q4-n 



i(27r)' 



q2 — -|- ie 



{q — — m2 + ie 

3r 



z(27r)" 



— ml + ie 



B2i{p,mi,m2) + g^iuB22{p,mi,m2) 



{q — p)"^ — m'2 + ie 



i {2Try 



g2 — 77^2 _|_ 



(q — pY — m'2 + ie 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



where Q is the renormalization scale and we regularize by integrating in n = 4 — 2e dimensions. 
The expression for ^0 can be integrated to give 



ylo(m) 



1 



m 



(B.5) 



where l/e = 1/e — 7£;-|-ln47r. 

Our A and B functions differ from those of Ref. since we use the Minkowski metric. Also, Ao, Bi and B22 differ 
by a sign. Equations (B.2 - B.4) contain an abuse of notation. The first argument of B-functions is the square root of the 
scalar p^, whereas elsewhere the p represents the external momentum four-vector. 
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The function Bq can be written in the form 

Bo{p,mi,m2) = T - / 
e Jo 



^ (1 — x) ml + X 777,2 ~ ^(1 — x) — ie 

dx In 



It has the analytic expression 

Bq{p, mi, 7712) 

where 



In 



P2 



fB{x+) - fB{X-) , 



X± 

^2 ^2 



s ± \J — 4p2 (777,2 _ 



/b(x) = ln(l-a;)-xln(l-a;-^)-l , 



and s = — 7/7-2 + ^1 • 

All the other functions can be written in terms of Aq and Bq. For example, 



7771, 7772) = 



1 

2^ 



-4o(7772) - ^o('77l) + (jP' + m\ - 777|)So(p, 777i, 7772) 



and 



522(f', 7771, 7772) = ^ I ^^^ol^l-l) +-4o(77l2)^ + (^777i + 777^ - ^P^^ Bo(P,V7ll, 7772) 



+ 



2 2 



^o(?"'2) - ^0(^^71) - (777! - ml)Bo{p, 7771, "72) 



1 



I 2,2-^2 
+ 777i + 7772 - -p 



We also define 

F{p, 7771,7772) 
G{p, 7771,7772) 
H(j), 7771,7772) 



^0('T7l) - 2^0(7772) - (2p2 + 2777i - 7772)-Bo(p, 7771, 7772) , 

{p^ - ml- ml)Bo{p, 7771, 7772) - Ao(mi) - Ao{m2) , 
4B22{p,mi,m2) + G{p,mi,m2) , 



B22(p, 7771, 7772) = B22ip,mi,m2) - ^Ao(777i) - ^^0(7772) . 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



(B.IO) 

(B.ll) 
(B.12) 
(B.13) 
(B.14) 



The functions F and G arise in scalar self-energies, with either a vector boson and a scalar or fermions 
in the loop, while H and ^22 occur in vector-boson self-energies, with either fermions or scalars in the 
loop. 

Appendix C: The gauge couplings 

In the remaining appendices we denote = sin^^vK; where 9w is the DR weak mixing angle, and 
= sin^^H^ = 1 — M^/M|, where 9w is the "on-shell" weak mixing angle and Mw, Mz are the 
gauge-boson pole masses. 
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The DR electromagnetic coupling is given by 



a 



l-Aa 



137.036 



wliere[^ 



Ad 



0.0682 ± 0.0007 

2 



Of 



27r 



\MzJ VM^y VM^y 3^ 



In 



(C.l) 



(C.2) 



and J2u indicates a sum over u,c,t, and similarly for J2dJ X^e- this expression, the number 0.0682 
includes the two-loop QED and QCD corrections given in Ref. ||4l|, as well as the five-flavor contribution 
AaJ2d(^l) = 0-0280 ± 0.0007 of Ref. ||. 
The DR weak mixing angle is given by [43| 
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where p is defined to be and JvB denotes the nonuniversal vertex and box diagram corrections 

given below. The W and Z gauge-boson self-energies are given in Eqs. ( D.4| ) and ( p.9| ). We compute p 
via m 
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pMl 
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We deduce the leading two- loop standard model corrections to Af and Ap from Ref. |41], 
Ar 



2— loop 
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2— loop 



47rs^c2 TT 
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where xt = SG^m'f /8tt'^\/2 and is the standard model Higgs-boson mass. For r < 1.9, p^'^\r) is well 



approximated by 
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^The coefficients of \n(Mw / Mz) in the expressions for Aq in Refs. ^ are both incorrect. 



31 



while, for r > 1.9, we use 
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For the case of the MSSM, we replace the function p^'^\mip/mt) with 
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We have not computed the corresponding Gj^mf higher-order contributions from the heavy Higgs bosons, 
but we know that they must decouple. Using the ansatz 
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Heavy Higgs 



/ sin a 
\sin/? 



we find these contributions are negligible. We do not include them in our results. 

The nonuniversal contribution to Ar is made up of two parts, one from the standard model and the 
other from supersymmetry. 
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The standard- model part is given by the well known formula [43| 
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The super symmetric part appears in Ref. [^, and more recently in Ref. |11|. We include it here for 
completeness. It includes box diagram contributions, vertex corrections, and external wave-function 
renormalizations. We neglect the mixing between different generations of sleptons, and we ignore the 
left-right slepton mixing in the first two generations, in which case the right-handed sleptons cr, Jir do 
not contribute. We find 
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The wave-function and vertex corrections are 
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The corrections (5Z^, and 6Vf^ are obtained from these expressions by replacing e — > The x- 

fermion-sfermion couplings a- ct and 6 - r ? are listed in Eqs. (P^20m^), while the char gino-neutralino- 

XiJ Jj XiJJj ^ — ' ' ' 

W couplings a-o^+yir and b-o^+yi^ are defined in Eqs. (p.l2 - D.13 ). In these expressions, the Bq, Bi, and 
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Cn functions are evaluated at zero momentum 
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where Q is the renormalization scale. 
The box diagram contributions are 
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where the functions Dq and D27 are 
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We checked our box diagram calculation against Refs. |5|, |TT|, and we checked our formulas for bZ and 
bv with those of Ref. Here (there) the formulas are written in terms of the couplings corresponding 
to vertices with incoming (outgoing) charginos and neutralinos. To compare we must make the transfor- 
mation o^jj ^^jj> except for couplings involving the chargino and down-type fermions, which remain 
unchanged. Also, their x^^W coupling differs from ours by a sign. 

The effective weak mixing angle is given in terms of the DR weak mixing angle, s^, via 



sin^e^r 



Tie kp 



where [45] 



with 



and 



1 + 



V,{Ml) 



c Hz^ (M|) -nz^(O) 



+ 



ac 



TTS^ 



In c 



a 



47rs2 

l-6s^ + 8s^ 
4^2 



y,(M 
/(I 



I), 



nefix) 




2Li2 



1 



vr 



1 + X 

J-) 

2x I X 



— +\n\l + x) 



1 + — I - ( tan"^ y\ . 



(C.23) 
(C.24) 
(C.25) 



(C.26) 



Here y = \/x/(4 — x), Li2 is the Spence function, and liz^y is listed in Eq. ( p. 15 ). We do not include 

qlept 



here the nonuniversal vertex super symmetric contribution to sin^^^'??*. The largest contributions can 



be obtained from Ref. |46|. 



Appendix D: One- loop self-energies 



In this appendix we list all the relevant self-energy functions which allow us to determine the one-loop 
fermion, gauge-boson, and superpartner masses. We explicitly include all of the necessary couplings. 
We perform our calculations in the 't Hooft-Feynman gauge, in which the Goldstone bosons and the 
ghosts have the same masses as the corresponding gauge-bosons. The gauge couplings g' , g, and gs, and 
the Yukawa couplings Xf are all DR couplings. The neutralino mixing matrix N, the chargino mixing 
matrices U and V, the Higgs mixing angles a and /3, and the sfermion mixing angles 9f are described in 
Appendix A, as are the normalizations of the Yukawa couplings Xf. The self-energies are given in terms of 



the Passarino-Veltman functions Ao, Bq, Bi, F, G, H, and B22 listed in Appendix B, Eqs. ( |B.5 ^ B.14 ). 

To streamline notation we do not write explicitly the external momentum dependence of these func- 
tions, e.g., we write BQ{p,mi,m2) as -Bo("ii, 717-2). Throughout this appendix we write s for sin, c for cos 
and t for tan, so that s/3 = sin/3, C2et = cos26t, etc., and for the sfermion mixing angles, Cu = cos0„, etc. 
Sub- or superscripts / denote a quark or lepton, and q denotes a quark. Inside a summation the 
subscript or superscript u denotes all up-type (s)fermions, u, c, t, u,,, u^, u^, and similarly inside a summa- 
tion the script d denotes all down-type (s)fermions, d,s,b,e,iJ,, and r. The sum J2f^/f^ denotes a 
summation over (s)quark and (s)lepton doublets, and the sum Eg denotes a sum over (s)quarks. Some 
terms are zero, for example Xi, = 0, and terms involving the right-handed sneutrino are absent. 
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In the self-energies listed below the 1/e poles are canceled by counterterms which relate the bare mass 
to the running mass. So, in the following DR self-energies we implicitly subtract the 1/e poles. 

Z and W bosons 



The full one-loop MSSM gauge-boson self-energies appear in Ref. [|| and subsequently in Ref. || 
The supersymmetric contributions are listed in Refs. |2l|, 

The self-energies of the gauge-bosons can be separated into transverse and longitudinal pieces, e.g. 



p2 



(D.l) 



The physical gauge-boson masses are the poles of the corresponding propagators, which involve only the 
transverse part of the gauge-boson self-energy. 



M| = M|(Q) - 7^e^i^(M|) , 



(D.2) 
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Here Mz{Q) and Mw{Q) denote the DR running masses which are related to the DR gauge couplings 
and vev's, as in Eq. ( |A.2D . The gauge-boson self-energies are evaluated at the renormalization scale Q. 
The transverse part of the Z-boson self-energy is 
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Co/3 
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where the summation X)/ is over all quarks and leptons, and the color factor NI is 3 for (s)quarks and 
1 for (s)leptons. The notation Sap denotes sin(a — and Cafj refers to cos(a — 

The sfermion-sfermion-Z couplings can be written in terms of the weak neutral-current couplings 
defined in Eq. ( [A. 7 ): 



vju = 9fLCf-9fRs}, Vf22 = gj^cj-gf^s}, Vfi2 = Vf2i = {g^ + 9fR)cfSf . (D.5) 
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The neutralino-neutralino-Z-boson couplings are defined by 



(D.6) 



and analogous definitions hold for f^^ and gf^^. We write the Feynman rule for the xX^/j. vertex, where 
X is a chargino or neutralino, as —ij^iiapL + bpR), where Vl,r are the usual chiral projectors (1 ^75)72. 



aj.oJ.oz and 67+r 



2c 



The couplings involving the unrotated V'° and ^+ fields satisfy fe^o^p^ 
The nonzero a-type couplings are 

For an incoming x^ and incoming x^ we have 
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(Here and in the following formulae which specify rotations, we adopt the summation convention for 
repeated indices.) 

For the transverse part of the T^-boson self-energy, we find 
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where the summation J2 fu/fd ^^^^ quark and lepton doublets, and 
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Wf2l 



CuCd 



Wfi2 
Wf22 
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The neutralino-chargino-VF-boson couplings are 

fijw = l«x°x+i^l' + l^x°x+H/l' ' 9ijW = 2ne(blo^+ 



(D.IO) 
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We write the Feynman rule for the neutralino-chargino-VF^ vertex as —i^^^aVi + bVn), and the nonzero 
couplings are 

For an incoming Xi we have the couplings to mass eigenstates, 



9 

V2 



while for an incoming Xj we have the couplings 



Nik U*i b^oj+w , 



eg 



Finally, we write the mixed Z — j self-energy as 

= {Us'^ -W)B22{Mw,Mw) - 2{M^ + 2c^p^)Bo{Mw,Mw) 
+ Xl^ce/(5/L 4^22(?n/,m/) +/So(m/,m/) 
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Quarks and leptons 

The fermion masses are defined as the poles of the corresponding fermion propagators. They are 
related to the DR masses, rhf, by the self-energies, Tjf{p'^), as follows 



nif = rhf{Q) — TZeT.f{m 



(D.16) 



The DR fermion mass is related to the DR Yukawa coupling and vev as shown in Eq. ( [A.l| ). Care 
must be taken in evaluating the DR vev. After evaluating the DR gauge couplings g' and g as outlined 
in Appendix C, we determine the DR vev via 



v\Q) = 4 



M| + 7^e^i^(M|) 
9'^{Q)+9HQ) ' 
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where Q is the renormalization scale (the argument of the Z self-energy is the external momentum; it 
implicitly depends on the scale Q as well). 
For the top quark, Ilt(p^) is 
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The neutral current couplings gf are defined in Eq. ( A.7 ). 

We write the Feynman rules for the Xiffj couplings as —i{aVL + bV^) (for vertices involving the 
chargino and down-type fermions the Feynman rule is iC~^{aVL+bVR), where C is the charge-conjugation 
matrix). We define 



iffj 
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In the unrotated ■0'') basis, we have 
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where the quantum numbers Yy and are listed in the table of Eq. ( |A.8| ) . These couplings correspond 
to vertices with incoming neutralinos and incoming charginos. To obtain the couplings to the mass 
eigenstates Xi a-^id xf , we specify the rotations 



^xtff' 
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The couplings to the sfermion mass eigenstates are found by rotating these couplings (both a- and 6-type) 
by the sfermion mixing matrix, 



xff[ 

\ff'2 



'xff'L 

xf fji 



(D.23) 
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The self-energies Sj(p^) for the other up- type quarks and leptons can be obtained from the previous 
formulae by obvious substitutions. For the bottom quark (and similarly for all down-type fermions), one 
interchanges t <^ b, Ca ^ Sa, and C/j ^ S/j. 

Charginos and neutralinos 

The complete one- loop self-energies for charginos and neutralinos are given in ^ ; we present them 
here in a matrix formulation. For the Higgs-boson contributions, refers to H, h, G^, and A, while 
represents and H^. The G^ and G^ are the Goldstone bosons; in the 't Hooft-Feynman gauge 
their masses are equal to Mz and M^y, respectively. 

We now describe the full one-loop neutralino and chargino mass matrices, from which we determine 
the one- loop masses. The one- loop neutralino mass matrix has the form 

+ U^^Ap')+^^Io{p')) , (D.24) 



where 



Here M^o is the tree-level neutralino mass matrix of Eq. ( A. 3 ), and the S^^'^ s(P ) matrix corrections. 
They allow us to determine the one- loop masses and mixing angles for arbitrary tree-level parameters. 
The one-loop chargino mass matrix is as follows, 

- _ - , (D.26) 



where is the tree- level chargino mass matrix of Eq. ( |A.4| ). The elements of -Mj^^o and A^^+ contain 

DR parameters at the scale Q. In particular, they include corrections corresponding to replacing Mz with 



Mz, obtained from Eq. (D.2). Similarly, tan/3 in the tree-level matrices is tan l3{Q). The self-energies 
^L,R,S are also evaluated at the scale Q. 

To obtain the mass for a given neutralino or chargino, for example Xij we first evaluate the matrix 



of Eq. (p. 24) with the momenta = m?o- We then solve for the eigenvalues of that matrix. So, in 

Xi 



determining four neutralino and two chargino masses, we construct a total of six different matrices. 

We compute the mass matrix corrections by evaluating two-point diagrams with unrotated neutralinos 
or charginos on external legs, and mass eigenstates inside the loop. We obtain the couplings associated 
with these diagrams by the following method. The neutralino mass corrections involve the couplings 
aro... which we obtain from the various couplings OyO... (for incoming x^) by leaving off one factor of N*j^. 
The neutralino mass corrections also involve the couplings h 70 which we obtain from the couplings b^o 

Vk'" Xi 

(for incoming Xi) by leaving off one rotation Ni^. We obtain the couplings aj;+... which appear in the 

^k 

chargino mass corrections from the couplings a-+ (for incoming xf) by leaving off one factor of V*j^, and 

Xi 

we determine the couplings bj+ from the couplings &^+... (for incoming xf) by leaving off one factor of 

^k Xi 

Uik- 

For the neutralinos, we have the one-loop correction 
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S^j is obtained from E^^ by replacing the couplings a^o... with b^o...- The correction is given by 
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The chargino mass corrections are given by similar formulae, 
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A;=l n=l 

^r{p^) is obtained from E^(p^) by substituting aT+ with is given by the following 
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formula, 
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In these expressions, the color factor is 3 for (s)quarks, and 1 for (s)leptons. The ipf f couplings 
are listed in Eqs. ( D.20| - p.23 ), and the tpx^ i^'xW couplings are given in Eqs. ( p.7 -D.8, D.12| - p!T4 ). 
We determine the '>p~^x^'~^ couplings from the following equations, which apply for incoming xt ^ 



X-i Xj 7 



e 6; 



Xi Xj 7 



e 8, 
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where we write the chargino-chargino-photon Feynman rule as —i'j^iaVL + ^Vr). 

We next list the xX-Higgs-boson couplings. We write these couplings in the unrotated Higgs basis 
{si, S2), {pi, P2), and {hf , /i^). These fields are rotated to obtain the mass eigenstate fields. The {H, h) 
rotation is given in Eq. ( |A.14D , while for (G*^, A) and (G^, H'^) we have 
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We write the Feynman rules for the X x^k couplings as —i{aVL + bVu) and for x x Pk as (aVi + VPr) 



These couplings are symmetric under i ^ j and satisfy bjnjn^ 



nonvanishing a-couplings are 
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The couplings to incoming neutralino mass eigenstates Xi are 
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and likewise for pn couplings. The couplings to Higgs-boson mass eigenstates are found by rotating these 
couplings, 
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and likewise for the 6-couplings. 

We write the Feynman rules for the x'''X'''"ii6utral-Higgs couplings as —i^aVi + bVn) for couplings 
with CP-even s-fields, and {uVl + ^Vr) for couplings with CP-odd p-fields. These couplings satisfy 
Hti'ts^ = ''i'ti^t^r. a'^d Ht^tPr. = ^he uouzcro a-couphngs are 

The couplings to incoming xf are obtained from these as follows, 

'xtxU^ = ^^kU*ia^+^+^^, b-+^p^ = UikVjib^+^+^^ , (D.38) 



and the same rotations apply for the p„-couplings. To find the couplings to Higgs-boson mass eigenstates, 
we rotate these couplings by the angle a or /3, just as for the x^x^s and x°xV couplings in Eq. ( p.36| ). 

The x''x'''~charged-Higgs-boson vertex Feynman rules are written —i{aT'L+bT'R), where, for incoming 
tp^, we have 

g' g 

To obtain the couplings to chargino and neutralino mass eigenstates with an incoming neutralino x^, we 
rotate these couplings as 

°x°x^^ = ^ikU*ia^oj+,,+ , b^^^+ = N,kVjib^oj+j^+ , (D.40) 
while for an incoming chargino x.j'^ we rotate them as 

«x°x^i = ^ikV;ib^oj+j,+ , b-o^+^+ = NikUjia^oj+f,+ . (D.41) 
To find the couplings to charged-Higgs mass eigenstates, we rotate both a- and 6-couplings by the angle 
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Gluino 

The gluino self-energy appears in Refs. I, I, g 0. The physical gluino mass satisfies 

m-g = MsiQ) - ne^gimj) , (DAS) 
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where Q is the renormahzation scale. 



Squarks and sleptons 



We find the sfermion masses by taking the real part of the poles of the propagator matrix 



Det 



- M){pi) 
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The matrix formalism allows us to determine the one-loop masses and mixing angles for arbitrary tree-level 
parameters. In this expression, the M? - , = are the DR tree-level mass matrix entries given 



in Eqs. (A. 5, A. 6): all the entries contain running DR parameters at a common scale Q. In particular, 
the DR tree-level matrix contains corrections from the replacements M\ = + TleIl'^^{M^) 

and ruf ^ TJif = ijif +7^eSj(mj). (The arguments of these self-energy functions are external momenta, 
not the scale Q.) The H? ? , (i, j = L, R) are the sfermion self-energy functions evaluated at the scale Q. 

hi 3 

Of course, for the first two generations of sfermions, both the tree-level and one-loop contributions to the 
off-diagonal elements of the mass matrices are negligible. Note 11 ? ? ^ W- - because of the absorptive 

JRIL JLJR 

part, which contributes to the mass-squared at 0{a^). 
For ati squark we have 
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and similarly for a tR squark, 
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The off-diagonal self-energy is 
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Inside the sum J2fj the sub- or superscript / refers to (s)quarks and (s)leptons, and in the sum J2f„j the 



sub- or superscript u refers to up-type (s)quarks and (s)leptons. The gj are defined in Eq. (|A.7|) . The 
electric charges e/, hypercharges Yf and third component of isospin Jg are given in Eq. ([A. 81) . These 
results are equivalent to those of Ref. [^] in the limit g, g', and — > 0. 
For the xff couplings, we have defined 



with analogous definitions for the LL and RR couplings. The xf f couplings are listed in Eqs. ( p.20[ - 
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The Higgs bosons refer to H, h, G^, and A, and refer to H^, . The i7„-//„-sfermion- 
sfermion couplings involve Gn and -Dn/j ^-nd are given in the following table, 
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We write the Feynman rules associated with the CP-even-Higgs-sfermion-sfermion vertices as 
and list the couplings ^gjj in the following table. 



We find the couplings in the /i^2 sfermion basis via 





Sl 


S2 


ulul 
urur 

ULUR 


"^^^ a ca 

c yuR<^i3 

Ay. 1, 


-^^QulSp + V^Xurriu 
^ 9urSi3 + \/2A„m„ 


d-Rd-R 
did-R 


amg^^cp + ^Xj^md 
^ gdiiCf3 + V^Xdrud 


nMz 

nMz 
c 9dRSp 



SnfLfh ^SnfhfR 



SnfRfh ^SnfRfR 



Sf 
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(D.54) 
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we obtain the couplings in the mixed /l,/?/i,2 basis by omitting the left-most matrix on the right hand 
side of the above equation. The couplings to the CP-even Higgs-boson eigenstates {H, h) are obtained 
from the couplings to (si, S2) using the rotation 




S2fifj 



The couplings '^QOf f. and ^aj.j. vanish for i = j, while for i ^ j they satisfy ^aj f 
write the Feynman rules for these couplings for incoming fi as A. They are 



(D.55) 





GO 


A 




ULUR 

didR 


-^ifiSp + AdCp) 




AuCp) 
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We obtain these couplings in the /l,_r/i,2 basis by a rotation as described after Eq. ( D.54 ). 

We also write the Feynman rules for the charged-Higgs-sfermion-sfermion vertices in the form —iX. 
The couplings Xq+jj, and X^j+^ p are 



urcIr 
uidR 
urcIl 



H+ 



-^C2fi - XuniuSp + XdrudCp 


Ad il^sp + AdCp) 



These couplings are obtained in the /i^2 basis via 




^^S2/3 - Kruucp - XdrudSp 
-XumdC/s - XdmuS/s 
Ad ifJ'Cp - AdSp) 
K ilJ'Sf} - A 



H+ULdL ^H+ULdR 
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H+URdL ^H+undn 




(D.58) 



We obtain the mixed sfermion basis couplings to UL,Rdi,2 (^'i,2C^l,_r) by leaving off the left-most (right- 
most) matrix on the right hand side of the above equation. 

The expressions for H^;^. are obtained from I^i-t. by interchanging the indices t <^ b, replacing u ^ d, 
and substituting <-> s^. The self-energy of a charged slepton (sneutrino) is given by a formula similar to 
that for a 6-squark (t-squark), with the SU (3) correction set to zero and with the appropriate SU (2) x U{1) 
quantum-number substitutions. 

Higgs bosons 



The full one-loop MSSM Higgs-boson self-energies appear in Refs. 
masses are the subject of Refs. [ 



] . Corrections to the Higgs boson 



36]. We discuss the relations between the self-energies and the pole 



masses of the Higgs bosons in Appendix E. Here we list the self-energies. 
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The Higgs-boson contributions to the Higgs-boson self-energies involve the trilinear and quartic cou- 

phngs, which we denote A^o^/o^^^, ^h+h^s^^ ^nd XHnH„HmH^, ^h„h„h+h-' ^here the refer to the 
H, h, G^, and A Higgs bosons, and the refer to the and Higgscs. The G^ and G^ are the 
neutral and charged Goldstone bosons, which in the 't Hooft-Feynman gauge have masses Mz and Mw, 
respectively. 

For the two CP-even Higgs bosons, we have 



+ 



fd 



+ E 



Nil 
2c2 



{p^ - 4ml)Bo{md,md) - 2Ao(m<i) -|- Ao{m^J + AQ{m^^) 
9fL (c/^o (^/i ) + s/^o (m )^ + gf^ (^s^Aq (mj^ ) + c^^o (mjj^ 



t1 



2F{mH+,Mw) + 



F(mA,Mz) 



C2 



+ 4 



2F(Mh/,Mw^) + 



F{Mz,Mz) 



C2 



, '22 



2M^Bo{Mw,Mw) + 



MlBo{Mz,Mz) 



+ 9' 



n=l '- m=l 
2 r 2 

n=l m=l 



+ E 
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(p^ - 4m^)So(m„,m„) - 2Ao(m„) + Ao(m.«i) + ylolwiij 



2F{mH+,Mw) + 



F(r».i,i\/z) 
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2F{Mw,Mw) + 



F{Mz,Mz) 



I '22 

+ 



2M^BoiMw,Mw) + 



1 r 

n=l m=l 

2 p 2 



n=l m=l 
-, 4 



E 



E 



«J«22 Xi X,' 
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16^2 Us,s,{p^) 



J^^s^c^J 2F(Mh.,Mh/) - 2F{mH+,Mw) 



+ 



F{Mz,Mz)-F{mA,Mz) 



+ 7 



2M^Bo{Mw,Mw) + 



MlBoiMz,Mz) 



n=l m=l 
2 



E 



X,- 
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where is the number of colors, which is 3 if / is a (s)quark and 1 if / is a (s)lepton. The neutral 
current couplings Qf are defined in Eq. ( A.?] ). 

The XiXi-Higgs couplings /^^^^, gf^^^^, f^j,^^, and g^j,^^ are defined by 



Ski 



^XiXjSk ^XiXjSi + ^XiXjSk ^XiXjSi 1 



(D.62) 



and the a-y^iXjSk ^'^'^ ^XiXjSk couplings are defined in Eqs. ( D.33 ), ( D.35| -38). The Higgs-squark-squark 
couplings ^Hf fk ^hf fk ^'■^ given in Eqs. ( D.53| - D.54 ). 
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We write the Feynman rules for the relevant quartic Higgs couplings as —iX, and define Xu^HnHmHrn = 
H„HmHm- -^is^ necessary XH^HnHmHm couplings in the following two tables, 





SiSi 


S2S2 


S1S2 


HH 




Qs2 _ 2 


-S2a 


hh 


_ ^2 




S2a 


qOqO 


C2/3 


-C2/3 





AA 


-C2/3 


C2/3 





G+G- 


+ S^C2p 


£2 - 5^02/3 


-C^S2/3 


H+H- 


- S^C2f3 


+ S^C2/3 


C^S2/3 





AA 


H+H- 




■^s^ - 1 






"^''2/3 


''2/3 








H+H- 


'-2/3 


^^-2/3 
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For the couplings involving (si, S2), we obtain the corresponding couplings in the {H, h) eigenstate basis 
by the following rotations 



inHnH H ^HnHnHh 

^H„H„Hh Xh 



Cq, Sqi 

Sa C-a 



Xh Xh 
Xh Xh 
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We write the Feynman rules for the trilinear Higgs-boson couplings as —iX, and define Xn^HmSi 
gMz/i2c)XH„H,nsr We list the Xn^HmSi in the following two tables, 







HH 






hh Hh 




CfsiScl 


-sl) 


- Si3S2a 


C/3(3s2 


- 4) + Si3S2a -'2Ci3S2a " Sf3C2a 


S2 




-cl) 


- CpS2a 


■8/3(34 


- sl) + Ci3S2a '2S/3S2a " C/3C2a 




G'^G^ 




AA 




G+G- H+H- G+H- 


Sl 


C2f3C/3 




-C2f3Ci3 


-•S2/3C/3 


C213C13 -C213C13 + '^C^C/S -S2i3Cf3 + C^SfS 


S2 


-C213S13 




C213S/3 


•S2/3S/3 


-C213S13 C2i3Sf3 + 2c^Si3 S213S/3 - C^C/3 



(D.66) 



(D.67) 



To obtain the couplings involving (si, S2) in the {H, h) eigenstate basis, we rotate the (si, S2) couplings 
by the angle a, as in Eq. ( p.55| ). 
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The CP-odd Higgs boson A and charged Higgs H'^ self-energies are 



+ 



+ 



4 



fu 



fu ^ ^ 



clAo{mui) + slAo{mu^) 



slAoiniu^) + clAo{mu2) 

c/3 ^ s/3 y J 



9_ 
4 

1 



2 r 2 2 

2F{mH+,Mw) + ^-^F{mH,Mz) + ^F{mh,Mz) 



4 r 4 



2 



BoiMw,mH+) + Y.^AAH+H-MmH+) 



n=l 



1 



+ E 



ftjAG{m-+,m-+) - 2g±^m-+m-+Bo{m-+,m. 
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ma) - 2XuXdmumdS2i3Bo{mu,md) 

fu/fd 

2 

+ E E^/^l,+^,d,^o(m^„m^P 

fu/fd » J = l 



+ E^/ Ai4-^i3"fe,+ 



-C2/3 



fu 



i2„2 9 TU„u^ 



+ ^ 



+ E^/(A^4-V3 

slf^F{mH,Mw) + clpF{mh,Mw) + F{mA,Mw) 

+ -^F{mH+,Mz) 
c J 

+ e''F{mH+,0) + 2g^Ao{Mw) + -^Ao{Mz) 



+ 



Sfj ^ CIS 
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n,m=l 



+ 



2 T\/r2 

w 



Bo{Mw,mA) + E ^H+H-HOHO Ao{mHo] 



n=l 



n=l 
2 4 



+ EE 

i=ij=i 



fijH+G(m^+,m^o) - 2gijH+m-+m^oBo{m-+,m^o) 



(D.69) 



where (sa/s) denotes cos(a — /?) (sm(Q — /?)). The ^j^, gj^ are defined in Eq. ( [A.7| ), and the Ig are 



listed in the table of Eq. ( A. 8 ). denotes the number of colors, which is 3 for a (s)quark. 
The XiXjA couplings f^^, gfj^, /°.^, and 5°.^ are defined by 



fijA 





2 




(^XiXjA 


+ 


^XiXjA 



9ijA 



!7^e (bl 



XjA'^XiX]A] , 
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and similarly for the fijH+i 9ijH+ couplings. The a^.j^-A, ^XiXjA couplings are given in Eqs. ( p.34 - D.38 ), 
and those of the charged Higgs are listed in Eqs. ( p.39| - p.42 ). The Higgs-sfermion-sfermion couplings 
^Af fk ^H+f fk given in Eqs. (D.56- p^58D . The /l^r basis couplings ^^ff °^ ^'i' ( P-56| ) also 
apply in the /i_2 basis. 



Appendix E: One-loop Higgs boson masses 



In this appendix we will present the formalism necessary to obtain accurate Higgs-boson masses at 
the one-loop level. The tadpole diagrams play an important role in determining the masses. The one-loop 
tadpole contributions are listed in Refs. ^] . At any given order in perturbation theory, minimizing the 
scalar potential is equivalent to requiring that the tadpoles vanish. At tree level, we have Ti = T2 = 0, 
with the tadpoles given by the DR relations 



T 1 

vi 2 
T 1 

— = - -M|c2/3 + rnt + + B^icotp 

V2 2 

where the Higgs-sector soft supersymmetry-breaking potential is 



soft 



rnjj_^ \Hi\ + m 



H2 



+ 



(B^ieijH\Hi + h.c. 



(E.l) 
(E.2) 

(E.3) 



At one-loop level, the total (tree-level plus one-loop) tadpole must vanish, so Ti — ti = 0, T2 — 12 = 0, 
with 



16vr2 ^ 

Vl 



Y^2Nf\lA,{m,) + EE^/oM . 
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+ 



+ 



9 



Ni3{Ni2- Na tan 9w) 



2 Xi 



i=l 

35^ 



-Tie 



Aq (m~r 



2Ao{Mw) + 



Ao{Mz 



^2, 



+ i^(2A^{Mw) + Ao{Mz\ 



(E.4) 



and 



167r" 



t2_ 

V2 



+ 



8c2 
8c2 

4 



2 



^o("^A) + 2ylo(m'H+) ) + Y^o("i//+) 



iVi4(iVi2 -iViitan^H/) 



1=1 

3<7^ 



771- + 

2 xT 



-7^e 



2^0 (Mh/) + 



Vi2Uil 

Ao{Mz) 



Aoim~+ 



9^C2I3 

8c2 



2AoiMw)+AoiMz)] , 



(E.5) 



where is the number of colors, 3 if / is a (s)quark and 1 otherwise. The ^0 function is given in 
Eq. ( |B.5| ); c denotes cos^vK and cp = cos /3, etc. The matrices N, U, and V are described in Appendix 
A and the coupHngs ^sif f 1 ^S2f f given by Eqs. ( D.53 , D.54|) . 

The DR (tree-level) CP-odd Higgs mass is given by 7h\ = — i?/x(tan/3 + cot /?), and Eqs. (ET, |E.2|) 
allow us to solve for the DR parameter, /i^^ and the pole massf^, uia, 



tan2/?(^m|^2tan/?- ml-^ cot/3^ - M| - 'ReX^F^ziMz) 



— {m\^-m\\ - Ml - TleUlziMl) - TZeUAAim\) + 6^ , (E.6) 
C2/3 V / 



where ^n^^ = rnjj^ — ti/vi, TTTjj^ = Tnjj^ — t2/v2- The self-energies and 11^^^ are given in Eqs. (D.4) 
and ( p.68| ), respectively, and = sj^ti/vi + Cpt2/v2- 

Having determined the physical CP-odd Higgs-boson mass rriA, we are in a position to compute the 
remaining Higgs masses. The physical mass for the charged Higgs boson is 



= m\ + + Tie 



(E.7) 



^^In case mA is very close to Mz, there is an additional 0(a) correction to niA from the oflt-diagonal element of the 
CP-odd mass matrix. 
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where the VK-boson self-energy is given in Eq. (D.9), and the charged-Higgs self-energy in Eq. ( D.69| ). 



The CP-even Higgs-boson masses are obtained from the real part of the poles of the propagator 
matrix, 



Det 

where the matrix Ml{p'^) is 



, mf = 7^e(pf ) , (E.8) 



Miip^) = \ \ . . (E.9) 



In this expression, and are the Z- and j4-boson DR masses (Mf = + 7?.en^^(M^), = 
-|- TleXiAA^jn^A) ~ '^'^^ self-energies n^.g^. are given in Eqs. ( |D.59| - p^ ). At one loop, the angle 



a diagonalizes the matrix Als(p^) for some choice of momentum p^; we choose = m\. 
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